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Abstract 

The Quantum Transfer Matrix method based on the Suzuki-Trotter formulation is ex- 
tended to dynamical problems. The auto-correlation functions of the Transverse Ising 
chain are derived by this method. It is shown that the Trotter-directional correlation 
function is interpreted as a Matsubara's temperature Green function and that the auto- 
correlation function is given by analytic continuation of the Green function. 

We propose the Trotter-directional correlation function is a new measure of the 
quantum fluctuation and show how it works well as a demonstration. 
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1. Introduction 

The Quantum Transfer Matrix(QTM) method is a so powerful tool to study thermo- 
dynamic properties of low-dimensional quantum systems. It was proposed by Suzuki 
[ 1 ] and was demonstrated by Suzuki and the present author |01 for the XY quantum 
spi n chains. It was widely used to one-dimensional integrable quantum systems (j- 
Holl and to non-integrable systems ITlll . We extend this method to study a dynamical 
problem by using a Transverse Ising chain as an example in this paper. 

The Suzuki-Trotter (ST) transformation maps a cZ-dimensional quantum system to 



ad+ 1-dimensional classical system adding an extra Trotter dimension II 11 1 1 lM 1311 . The 
QTM T transfers along to the real dimension as shown in fig[T] The free energy of 
the relevant quantum system in the thermodynamic limit is given by the one maximum 
eigenvalue A max of T as follows 10, S . 

-j3f = lim \ogA max , where T\uj max ) m = A max \i[r max ) m . (1) 

m— >oo 

Where we add a suffix m to indicate the finite Trotter number m. The thermal average 



of a physical quantity Q is expressed as with the normalized eigenvector iff max IU4I1 



(Q) = lim (ifr max \Q\i// max ) m . (2) 
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Suzuki et. al.| 14] obtained a magnetization (tr r ) and static real-directional correlations 
(crfcTj) of the present model. 

The main idea to study auto-correlation functions is follows. Consider the Trotter- 
directional correlation functions of the ST-transformed classical Ising spins in the extra 
dimension 

(S oS r ) m - (l]/ max \S oS r\*l>max)m (3) 

where r is a Trotter-directional distance of two spins. This correlation function is also 
expressed with the original quantum spin cr as 

(S„S r ) m = Tr e- fm <^ r '> /m o*e- fmr " n o* /Z m 

— >(cr?(r)cr?) ( as m — » oo) (4) 

with 

of ( T ) = e^crfe-^, (a = x,y,z) and r = fir/m. (5) 

Thus the Trotter-directional correlation function can be interpreted as a Matsubara's 
temperature Green function 11511 after taking the limit of the Trotter number m — > oo. 
The auto-correlation function is obtained by analytic continuation of r to imaginary 
time it IU5I1 . 

We propose the Trotter-directional correlation function as a quantitative and con- 
crete measurement tool of quantum fluctuations. A quantum state is usually described 
as a super position of classical states such that a two spins singlet state is defined as a 
classical | fl) state minus | J,f) state. In the ST-transformed system we consider that the 
classical Ising states are stacked along to the extra dimension. The original quantum 
spin state is represented as a superposition of these stacked Ising spin states as shown 
in fig[T] When all Ising spins have the same states, i.e., (SqS r ) = 1, the quantum fluctu- 
ation is zero, otherwise the quantum fluctuation exists. If (SqS r ) = 0, the fluctuation is 
maximum. Thus the correlation function can be a measure of the quantum fluctuation. 

We study r (or r) and temperature dependencies of our correlation function in the 
§4 and we show it works well as a measure of the quantum fluctuation. We re-derive 
auto-correlation functions by the QTM method in §5 from the results in the §4 and 



compare with the known result 11161 1 1 81-422 ] 



2. QTM and the maximum eigenvalue 



The Hamiltonian of the present transverse Ising quantum chain ||23j |24J is defined 



(6) 



By the ST-transformation we have a two-dimensional Ising model shown in fig JT] 
Its partition function is ifll H [l4ll 



N m 



Z m = A*"Trexp £ J] -SijSmj + YJS ijS 



(7) 
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Figure 1 : The original quantum spin (cr) chain and the ST-transformed two-dimensional Ising (S ) system. 
The QTM 7~ transfers along to the real direction. 



where 



7 /l 2y\V2 

exp(-2y,„) = tanh — , A m = - sinh — 

m \l ml 



(8) 



The system length N is assumed to be infinite. The periodic boundary condition is 
required for the Trotter direction. We can see that the vertical (horizontal) interaction 
Y m {—) increases (decreases) as the temperature increases. 

We apply the exact solution of the two-dimensional Ising model obtained by Schultz 
et. al J25ll using the transfer matrix method 12611 to the present model. The QTM T is 
defined in a symmetric (and then hermitian) way as follows j^B- 



T = V, I/2 ViV2 /2 with 



I 2K\ m l 2 

V l = (2 sinh — J exp(^ J] S % ^ = ex P^'» Tj S PU^ 

j j 

exp(-2r,) = tanh-. (9) 
m 

Here the spin operators S X ,S Z are used. The original quantum spin cn is mapped to 
the Ising spin Sy in eq.(|7]i and is remapped to S x in their formulation 12511 . Then the 
correlation we want to evaluate is given by (SoS , ) = (SqS x ) and it is a z-component 
correlation function of the original quantum chain as given in eq.©. 
The maximum eigenvalue A max is given bv l25ll 



/ . , 2KV"' 2 ,l r v 
Amax = 1 2 sinh — I exp ( - ^ s 9 J, 

1 

2y u 2K 2y 2K 

coshe 9 = coth — coth cosec — cosec — cos q, 



mm mm 



3 



1 3 m - 1 

q = ±-,±-,---,± , (10) 

mm m 

with the help of Jordan- Wigner, Fourier transformation and Bogoliubov diagonaliza- 
tion. The associated maximum eigenvector is a vacuum of fermions. 

We note that the critical line of the ST-transformed two-dimensional Ising model in 
eq.(|7]i is given by 

2K 2K 2y 

1 = sinh — sinh 2y' m = sinh — / sinh — (11) 
m ' m m 

as the ordinary two-dimensional Ising model. This means that our system is critical 
when y — K without the dependence of the Trotter number m i27ll . When y > K the 
system is disordered and when y < K it is ordered 02 ll 12411 . 



3. Trotter-directional correlation functions 



The correlation function (SoS r ) m at a finite temperature is expressed as an r x r 
Toeplitz determinant 123142511 . 



{SqS r) m — (lf/ m ax\S qS r \lf/ max )„ 



go 


g-l g-2 


gl-r 




gi 


go g-i ■ 


■ g2-r 




g2 


gl go 


g3-r 


(12) 


gr-\ 


gr-2 gr-3 ' 


■ go 





It becomes zz-correlation function of the original quantum Transverse Ising chain : 
{cr z {f)cr z ) = lim,„^ co (5()5 r ) m . The matrix element g^ is an expectation value of two 
fermions apart k + 1 position and it is denoted as an in the reference [25]. 

r 1 in —7(v4-fTi/m 1 



(—\) k -i.fi - e iq e~ 2iy+K ^ m l — e~'9e2(y-AT)/»n ^1/2 



(-If 



s 

q>0 



I — e -iq e -2(y+K)lm J _ £ iq e 2(y-K)/n 
_ „2y/m 



2 cosh(2K/m) cos kq - e^ /m cos(fe + l)q- e^ 1 " 1 cos(k - \)q 

TJ2 



[(cosh(2y/m)cosh(2/T/m) - cosg) 2 - (sinh(2y/m) sinh(2^T/m)) 2 j 



for K > y. With the help of an integral formula 



1 r 1 

n Jo a + b > 



g ± k becomes 
(-\) k 



g±k 



SI ■ 



■d6= , 
cos6» Va 2 - b 2 



(a > b), 



(13) 



(14) 



-de. 



2 cosh(2K/m) cos kq - e lylm cos(k +Y)q- e~ zr/m cos(Jt - l)q 
hut ^ Jo cosh(2y/ra) cosh(2A"/w) + sinh(2y/ra) sinh{2K/m) cos 6 - cos g 

(15) 
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The integrand can be simplified by using formulae of trigonometric functions and 
£o>o cos nq = 0(n + 0). With the following s and s' defined as 



2y 2K 2y 2K 

cosh s = cosh — cosh H sinh — sinh — cos 6, 

mm mm 

2y 2K 2y 2K 

s sinh s = sinh — cosh h cosh — sinh — cos 6, 

m m m m 



(16) 



the integrand becomes 



m sinh — | s' sinh ks + cosh ks ) + } — - sinh s{ s' cosh ks + sinh ks 

J cos q - cosh s 



m \ 

sinh(2y/m)/ 

- - s sinh( 
)V 



2 sinh(2y/m) 



cosh(mi/2)\ 2 
Here we have used the following formula 

1 



q>0 

k)s + cosh(— - k)s\. 



sis' . 



) 

(17) 



Z 

9 >o 



cos g - cosh ^ 



w tanh(ms/2) 
2 sinh i 



Thus we have the matrix elements 

. ,sinh(2y/m) f 1 
n Jo cosh(« 



m m 

s sinh( k)s + cosh( k).? 

ms/2)l 2 2 



2y sinh(2-y/ra) 

go = cosh 

m 7T 



, , ms , 
s tanh — d8. 
2 



(18) 

d0,(k > 1) 
(19) 



These are valid for the both cases of y > K and y < K. We note here that g% satisfies 



gk(y) = g-k(-y), 



g m . k (y) = (-iy»- V*(r) (k * 0). 



(20) 



We can easily check that (SoS m ) m = (SoSo)m = 1 is fulfilled. 

The magnetization (<x A ) of the original quantum spin chain is related to the nearest 
neighbor correlation function (SqS \) m = go as follows II 1 411 . 



2y 2y 
{a 31 ) = lim coth go/ sinh — = lim 

m— *oo in in m— >oc 



i r > u ms 

s tanh — 
Jo 2 



c/(9. 



For large m, s and s' in eq.dToTl becomes 
2 



Jy 2 + K 2 + 2yK cos 0, 



y + K cos 
yjy 2 +K 2 +2yK cos, o' 



(21) 



(22) 



and thus we reproduce the exact solution. 

The diagonal correlation function is expressed as a four-body Ising spins correlation 
as follows. Using a Ising spin state \S ) and its completeness £(si \S)(S | = 1, we obtain 



(<T x (T)a x ) = lim V (S\ 

m— >oo / J 



lim 



(S,S',5",S"'| 
1 

sinh 2 (2y/m) 



cosh — - 5 r+ i5 r )(cosh — — S \So) 

HI tn 1 



\S"){S"\e TrH \S'"){S'"\e 2y '»< s ' s " 



2y 



\S)/Z m 



(23) 
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Here we have used an identity, 

(S \& t e 7 °* ,m \S ') = e- 2Y »' SS ' (S \e 7 °* /m \S '). (24) 
The 4-body correlation of Ising spins on the same Trotter axis is expressed by gk as 

(S S \S r S r +\)m ~ 8q~ grg-i- (25) 



Thus we have 



<cr*(T)cr*> = lim ————( cosh 2 -1 - 2g\ cosh -2 + - grg _ r )\ 
sinh z (2y/m)\ m m I 

= (c^) 2 - lim —J-— grg _ r , (26) 
m->oo sinh (2y/m) 

Similarly, the yy-correlation function is given by 

{a y (T)a y ) = lim l(-iS r )( cosh — - 5 r+ i5,.)(-/5 )( cosh — - 5i5 

m-^oo s inh 2 (2y/m) \ v m y v m 



lim , 



- 1 r t 2y 2y / \ 

(cosh — + l)<5 5 r >„, - cosh — \{SoS r+ \) m + {S\S r ) m ) 

(27) 



sinh (2y/m) 



and which satisfies the following identity 1 28] for the present model at the limit m — > oo. 

^(o*(t)<j*) = r 2 <^(T)c^>. (28) 

4. Evaluation of the (SoS r ) m correlation functions 

We study the r and the temperature dependencies of the correlation functions {SoS r )n 
in this section. We cannot apply Szego's theorem ll29tl to our Toeplitz determinant in 
eq.(fT2l>. If we take m — > oo for a fixed r, all gk becomes except go — > 1 and thus we 
always have a wrong result {SqS r ) m - 1. We need to evaluate the determinant with 
a finite m and then take m — > oo limit at final to get a correct result of the original 
quantum chain |2]. 

Define 6 = r/m to indicate the normalized distance of two spins. We assume m 
and r are very large but finite. We take a perturbative approach for the small 6 in 
the first subsection. We do numerical computation for the general 5 in the successive 
subsections. 

4.1. Analytic approach for the small 6 

We calculate the determinant for the small 6 up to the second order analytically. 
Expand the matrix elements gk up to k 2 order, 



sinh 22 



m 



, 1 r > u ms 

1 s tanh — , 

n Jo 2 



de 

o 

(29) 



If*, mi £ 2 1 r , 
+Jk- j(i + tanh — ) c/0 + s (±1-5 tanh 

Jo 2 2 7T J;, 



ms 

— de 



6 



In general the determinant of a Toeplitz matrix whose elements has up to k 2 dependence 
can be obtained with the help of the standard matrix procedure and is given in Appendix 
A. 

After tedious but straightforward calculation we have 

(S S r ) m = 1 + ci(m)d + c 2 {m)5 2 + 0(6 3 ) 

with ci(m) = -2y/ - (2y 2 / () - lyh +h) + 0( m - 4 ), 



c 2 (m) = 2y 2 - ^4(/ 2 - hh + K 2 ) + 0(m- 4 ). (30) 



3m 

r 

3m 

Where integrals are defined as 

i r y+Kcoso n — : 

h = - ; tanh Jy 2 + K 2 + 2yK cos 6 dO = <cr x >, 

* Jo Vr 2 + K 2 +2yK cos 6 y 

h = ^ J ^y 2 + K 2 + 2y^ cos 6» tanh yjy 2 + K 2 +2yK cos d6, 

I 2 = - J^(y + K cos 6) yjy 2 + K 2 +2yK cos 6 tanh yjy 2 + K 2 + 2yK cos 6 d6, 

1 (yK 2 (y 2 + 2K 2 + 3y^cos 6») sin 2 ri " 

7 3 = - - — V , = tanh Jy 2 + K 2 + 2yKcos6 

n Jo V Vr 2 + ^ 2 + 2r^cos6' 3 v 

+ r 2 K 2 (y + Kcose)sin 2 e v < ^ 

(y 2 + /T 2 + 2y^ cos 6») cosh 2 ^y 2 + K 2 + 2yK cos9 ' 

The first integral 7o is nothing but the magnetization (cr A ) = --^(JSf) of the original 
quantum chain. It takes a finite value between and 1 depending on the ratio K/y in 
the low temperature region and it becomes y at the high temperature limit. 

We can check easily that the result eq.d30i> agrees with a perturbative calculation of 
the temperature Green function for a small t - /36 , 

(o- z (t)o- z ) = ((<r t + ['H,<r t ]r+^[ < H, [ < H,o z ]]t z + ---)cA = 1 - 2r<cr r )r + 2r 2 r 2 + ■ • • . 

(32) 

4.2. Numeric evaluation for K/y — 0.7 in the disordered region. 

To study the behavior of our correlation functions for the whole range of the pa- 
rameters 6 and the temperature we perform a direct numerical computation of the 
integrals of eq.([T9il and the determinant of eq.(fT2l. The data are taken for K/y - 
0.7, 1.0, 1.3, 2.0, m — 8 ~ 1024 and y = 0.1 ~ 64 while y/m < 1. The distance r varies 
r = 1, 2, 3, • ■ • , m/2, m/2 + 1. Due to the periodicity for the Trotter direction, our cor- 
relation function has the same value at r and at m - r. It is a good test for our numeric 
computation whether the function has the same value at r = m/2 - 1 and r = m/2 + 1 
or not. Our data has 9 digits accuracy even in the worst case. 

Figure |2] is the 6 dependence of the correlation function for the high temperature 
K = 0.7 and y = 1. We omit the data for 5 > 1/2. It takes the lowest value at 5 — 1 /2. 
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Figure 2: (SoS, ) m versus S = r/m at K = 0.7 and y = 1.0 for various m = 8 ~ 1024. The fitting curve for 
m = 1024 is also drawn. 



Its ra-dependence is very week as was shown in the analytic evaluation in eq.(f30b (m 2 
dependence) so that we can not distinguish each curves of the different m in this figure. 

We adopt the following two parameters (b and c) fitting curve to analyze the 6 
dependence, 

MS) = 1 -(cosh(c/2)-cosh(c( 1/2 -£)))/ cosh £ (33) 

which satisfies /#(0) = 1. This curve is exact for K — as shown in Appendix B. 

(S S r ) m = cosh(y(l - 2c5))/ cosh y. (c = 2y, b = y) (34) 

The fitting gives excellently good agreement to the numerical data as shown in figE] 

At a low temperature the correlation function shows exponential decay as shown in 
fig|3] We use a new fitting curve which does not satisfy fi(0) = 1 any more to estimate 
the main term. 

f L {5) = cosh(c(l/2 - 5)) I coshfo. (35) 

The results are 

cosh(21.2(l/2-c5))/coshll.8 for y = 32, 
cosh(40.6(l/2-c5))/cosh21.8 for y = 64, (36) 

and are drawn in fig. [3] The values of b and c for y — 64 are about twice larger than 
those for y = 32, respectively. 

To study the temperature dependence we plot the data for 6 — 1/2 in fig |4] The 
exponential decay for the y is clearly shown. By fitting the exponential function to the 
data of m = 32, 128, 512, we estimate the gradient approximately 

exp(-0.32 r ). (37) 
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Figure 3: Semi-log plot of correlation functions versus 6 for y = 2,4,8, 16,32,64 from top to bottom. 
K/y = 0.7 and m = 512. Two thin lines are fitting curves of eq.^36}. 




Figure 4: Semi-log plot of the correlation function at S = 1 /2 versus y, K/y = 0.7. m = 32, 128, 512. Lines 
are fitting curves shown in the legend. 
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Figure 5: Log-log plot versus S for various y. K = y and m = 1024. The line for y = 64 is a fitting curve 
given in eq. j38t . 



From these results we conclude that the correlation function decays exponentially 
to zero for both the 6 and the temperature (y) in the disordered region. 

4.3. Critical region: K — y 

We show two graphs figs. [5] and [6] for m = 1024. The former shows the 5- 
dependence and the latter shows the y-dependence. The both graphs clearly show a 
power law behavior of the correlation functions in the low temperature. The fitting line 
for y = 64 in fig [5] is 

1.01 xr 024 . (38) 
The fitting line shown in fig|6]for r — 512 (S — 1/2) is 

0.72 x y- (n5 . (39) 

These exponents remind us the two-dimensional Ising model which has rj — 1 /4 and 
£=1/8. 

4.4. Ordered region: K/y = 1.3 

The big difference from the other regions is that the correlation function decreases 
and saturates to some value as the temperature decreases in the ordered region. In fig JV] 
the saturated value is about 0.8 for K/y — 1.3. This means the quantum fluctuation is 
small in this region. The increase for the large 6 and the low temperature (y > 2) is due 
to the finite m. As is seen in figH] this increase disappears as m — > oo. 

We estimate the saturation values as follows. Plot the data of 6 — 1/2 versus 1 /m 
in figlH] Fit them with a quadratic equation of l/m and find the limit values of m — > oo 
for each temperature. Three these fitting curves are also drawn in fig|8] The three limit 
values almost coincide to each others such that 0.79934 for y = 32, 0.79928 for y = 48 
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<SoSr> K = y . m=1024 




0.1 1 10 y 



Figure 6: Log-log plot versus y for various position r. The line is a fitting curve for r = 512 given in eq. )39t . 




Figure 7: Normal plot of m = 512 for the several temperatures. 
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Figure 8: Saturation values versus m ' for K/y =1.3. y = 64 ~ 32 from top to bottom, m = 
224, 256, 320, 384, 448, 512, 640, 768, 896, 1024. Fitting curves are quadratic equations of 1/m. 



and 0.79912 for y — 64. Thus we conclude the saturation value for K/y - 1.3 at the 
zero temperature is about 0.799. Additionally, when K/y = 2.0 the saturation value is 
estimated about 0.930. 



These values are equal to the square of the order parameter (cr z ) 2 — (1 - y/K) l/li at 



T = |21|,|22|,|24H. 

When y = (K/y 
go = 1- 



oo), the value is (SqS r )m = 1 f° r an y r since g 



±k 



and 



5. Auto correlation functions 

We assume that the analytic continuation from the temperature Green function to 
the time-dependent auto-correlation function is achieved by just replacing t with the 
imaginary time it in the present case. This assumption is confirmed by deriving the 
same result as the known results lfl6l 18-22]. 



5.1. {(T Z {t)(T Z ) 

When t =! 0, replacing 6 — > t//3 — > it/fi in eq.d30b for m = oo, 

(trXOcr 2 ) = lim (S S r ) m = 1 - Ur(o*)t - 2T 2 t 2 + 0(t 3 ) (40) 

m— »oo 

This result agrees to Brandt and Jacoby[ 18] and Perk et.al il9ll with rescaling t by t/2 
due to the definition of the Hamiltonian. 

We adopt the following functions for the general t with many fitting parameters to 
fit whole the range < 8 < 1/2 of the numerical data shown in figs l3l5l andl7l Each 
fitting function is chosen to have the main exponential or power decay terms found in 
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Figure 9: Time dependence of the auto-correlation functions (a- z (t)cr z ). Solid lines are real part and dotted 
lines are imaginary part of the functions. The graph of y = iO,K = l should be compared to figs. 3 and 4 of 
Perk and Au-Yang [22], and the graph of y = 10, K = 13 to figs. 6 and 7 of it. 



§4. 5 is replaced to t = [55 (ft = y). The used data are for m = 512 and y = 10. The 
results are 

3 le-° MT 

0.028 + 0.48e" 2 ' 7T + 0.64e~ a74T + for y = 10, K = 7, (41) 

t - 20. 

- 2.9 + 8.3(1 - 1.0e 04lT )-r 1/4 + 0.30e~ 2 - 7T + 3.7e 048T for 7 = K = 10, (42) 

0.088e~ 2 - 2T 270e 0032T 4700e" 25T 

0.85 + T + r fory= 10,^= 13 (43) 

0.44 + t (70 + r) 2 (400 + r) 3 r ' 

and are shown in figs|9]after replacing t to it. These shows almost the same character- 
istic behavior obtained by many authors 1118 . Sin such that a fast+slow oscillating 



decay around zero for y > K, a simple oscillating decay for y — K and a saturation to a 
non-zero value with a oscillation for y < K. It is, as is well known, not easy to get ac- 
curate results by an analytic continuation from numerical data. We need to consider m 
dependence to get more accurate result as has done in §4 or need a more sophisticated 
method. It is instructive that the Trotter number m = 512 is large enough to estimate 
the leading singularity as shown in §4, however, it is small to estimate the whole range 
of the time dependence of the auto-correlations quantitatively. The present result is an 
example obtained by a simple extrapolation. 
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5.2. {a- x {t)a- x } 

Substitute g r and r = it = /3r/m into eq.(l26b and take the limit m — > oo we have 



, [i r n r + jcos0/ \ 

'— — I tanhySw cosh(2/w?) - sinh(2/w?) )d6 

[n Jo w \ I 

— J" ^cosh(2/wf) - tanhfJw sinh(2iwt)Jd8 



(44) 



for a finite temperature where we have abbreviated w — VF 2 + J 2 + 2TJ cos 6. Thus we 
have obtained the Niemeijer's exact result IKj 17, 2(J by extending the QTM method. 



6. Summary and Discussions 

We have studied the Trotter-directional correlation function which is interpreted as 
a Matsubara's temperature Green function by the Quantum Transfer Matrix method. 

By applying analytic continuation to the correlation function we examined the 
auto-correlation functions {cr-{t)cr z ) and (cr x (t)cr x ). The numerical result for the zz- 
correlation caught the characteristic behavior of the time dependence of them. We 
have re-derived the joc-correlation exactly. 

Our formulation are based on the ST-transformation and thus f3/m should be small 
enough to get a meaningful result of the relevant quantum systems. In practice we 
use a discrete t - fib - f3r/m to estimate the continuous r or t dependence of auto- 
correlations. Thus we need a large m to get more accurate results in numerical calcula- 
tion. 

We have demonstrated that the Trotter-directional correlation function is a useful 
tool to measure the quantum fluctuations. It has the same value of the square of the 
order parameter (cr z ) 2 at T = 0. The shrinkage of a spin length which is given as 
the order parameter is often used as a measure of the quantum fluctuation. While the 
order parameter (cr) is zero for T > for the present one-dimensional system and then 
it can not be used as a measure, our correlation has a finite value which varies from 
(<t z ) 2 (T — 0) to 1 (T — oo). When y — which means no quantumness (1 dimensional 
classical Ising model) , our correlation is equal to one and thus indicates no quantum 
fluctuation correctly. The square of the magnetization {cr ) 2 - (ctoctr) (R — > oo) is 
considered as a projection to an another spin in distance R. Our correlation (cr:(T)cn) 
is a projection of itself and also it shows how the fluctuation increases as a function 
of "time"(T). Thus our correlation function and the order parameter are complemental 
measures for the quantum fluctuation. 

We emphasize that the Trotter-directional correlation function can be defined for 
any ST-transformed classical systems mapped from quantum systems. What we should 
observe is the classical Ising spins {S } and thus we can easily done by the Quantum 
Monte Carlo simulation or others. By using many Ising spins which are on different 
stacked layers or different real positions, we can study general many temperature/time 
Green functions. 

Finally we mention Bethe-Ansatz systems. The maximum eigenvalue and the as- 
sociated eigenvector of the Heisenberg chain are well known yl 0, [l^l so that the 
application our method to them is a future problem. 
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Appendix A: Determinant of a Toeplitz Matrix 

The determinant of the rxr Toeplitz matrix with the elements 



[ a + (i- j)b + (i- j) 2 c 
M u = I X 

{ d + (j-i)e + (j-i) 2 f 



1 < j <i<r 
\<i=j<r 
1 <i< j <r 



(45) 



\c-f. 
(Ao) r - 1 



is expressed as follows. 

detM. (^jW-f^jW 

{ - 6c 2 fA QP (r + 1) - T iP {r) + T lP (r - l)}/2 
-rc|(zT + 4c/(c - /))/j(0 - zT>(r - 1)) 
+r 2 c 2 f(c-f)\p{r)-p(r-\)} 



(c-ff 



(Ao) 



■-2 



(C-/) 3 



6c/ 2 A^(r + 1) + A z 3 <?(r) - A z 2 <z(r - 1) f/2 



+^{Ao(zi - 4c/(c - /)) 9 (r) - A QZl q(r - 1)} 



+r 2 c/ 2 (c-/){A 0? (r)-Ao< ? (r-l) 



(46) 



Here 



zi = c 2 e - / 2 £ + 2£c/ + 2c/(c - /), 

z 2 = c 2 e 2 + 2c(-cd + (b + c)e)f + (b 2 + 4-bc + 2c(2a - 3X + c + 3c/))/ 2 - 2(a + b + c)f, 
zi=zi- 4/zi + 2c/ 2 (3a - 3X + 3£ + 3e + c - /), 

A =X -d, A\ = -3X + a + b + c + 2d + e-f. (47) 

The symbol overline z means an exchange (a, b, c) <-> (d, e, f) of z. (transpose of the 
matrix) And 



p(k) = 



1 



q(k) = 



(a-/3)(J3-y)(y-a) 
1 



(a-/3)(/3-y)(y-a) 

pTk) = (a /m) q(k). 



a k+l (J3-y) + /3 k+l (y - a) + y* +1 (a - j3) 
(a/3) k (a + (j3y) k (J3 -y) + (ay) k (y - a) 



(48) 



a,/?, y are roots of Aqx 3 + A\x 2 - A\x - Ao = . 

Multiply (-1)' x (-1) ;_1 = (-l)' -7-1 to Mij to suite our matrix M' with elements: 



' (-iy-J- l {a + (i-j)b + (i-j) 2 c} 



-Xa 



(-l)H-i{d + (j-i)e + (j-i) 2 f} 



1 < j < i < r 
\<i=j<r . 
1 <i < j < r 



(49) 
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The relation M and M' is 

detM' = (-l) r detM. 



(50) 



Appendix B: {SqS r ) m for K — (derivation of Eq.i34j) 

When K — 0, the system reduces to a one-dimensional Ising chain for the Trotter 
direction with a spin interaction y' m . The correlation function is simply given by 

tanh f y ;„+tanh'"- r y ;, cosh(y(l - 2£)) 

v>(K>/>m = 7— — r^— ; = r (51) 

1 + tanh y' m cosh y 

under the periodic boundary condition with using eq.©. 

Alternatively, we can derive this result from eq.dTZll. The s and s' of eq.dToTl are 
now 

s=— , s' = l. (52) 
m 



The elements gu are simplified as 

k _ x sinh(2y/w) 



(-if 



cosh y 



/ 2k \ 
+ expl +(1 ))y 



2y 2y 

go = cosh sinh — tanhy. (53) 

m m 

Thus the matrix of eq.(fT2b can be diagonalized and the determinant is obtained directly. 

As the third method, we can derive the result eq.d5"TTl without the ST-transformation 
by using a Heisenberq equation: 

—cr z (f) = [<H,o- z (T)}, cr(-r) = e' w V z e^ T (54) 
dr 

where r is imaginary time and H = — J^j Tcrj. Using twice this equation, cH^t) follows 
a differential equation and it becomes as 

d 2 

-T^o^(t) = rV'(r), ^(t) = c x e 1YT + c 2 e- 2FT . (55) 

ci and C2 are constant 2x2 matrices and are defined by an initial (o- z (Q)cr z ) = 1 and a 
periodicity (cH^rtcr 2 ) = {cr z (B - t)o~ z ) conditions! 15]- We have again 

cosh(y - 2Tt)) 

(^(t)^) = lim (S S r ) m = ^— ^. (Ft = yd). (56) 

m->oo cosh y 
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